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The recently observed accelerating cosmological expansion or "dark energy" may be either a 
negative pressure constituent in General Relativity (Dark Energy) or modified gravity (Dark 
Gravity) without any constituent Dark Energy. Low- or high-curvature modifications of Einstein 
gravity are distinguished by the spacetime (Ricci) curvature of their vacua. If constituent Dark 
Energy does not exist, so that our universe is now dominated by pressure-free matter, Einstein 
gravity must be modified at low- curvature, becoming asymptotically de Sitter. 
The dynamics of either kind of "dark energy" cannot be derived from the homogeneous expansion 
C ~ ) , history alone, but requires also observing the growth of inhomogeneities. Present and projected 

C _~) • observations are all consistent with a small fine-tuned cosmological constant, but also allow nearly 

£N| ' static Dark Energy or gravity modified at cosmological scales. The growth of cosmological fluctua- 

tions will potentially distinguish between static and "dynamic "dark energy". But, cosmologically 
distinguishing dynamic Dark Energy from Dark Gravity will require a weak lensing shear survey 
more ambitious than any now projected. Dvali-Gabadadze-Porrati low-curvature modifications of 
Einstein gravity may also be detected in refined observations in the solar system or in isolated 
£SJ ' galaxy clusters. 

We review local and cosmological tests of General Relativity and modified gravity. Dark Energy is 
epicyclic in character, requires fine-tuning to explain why its energy density is just now comparable 
^v^J , to ordinary matter density, and cannot be detected in the laboratory or solar system. This, along 

with braneworld theories, now motivate searching for Dark Gravity on solar system, galaxy cluster 
qq 1 and cosmological scales. 
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I. INTRODUCTION: COSMOLOGICAL SYMMETRY VS. DYNAMICS 

The greatest mystery in cosmology is the "dark energy" source of the late (rcdshift z < 1/2) cosmological accel- 
eration. This "dark energy" may be static or dynamic and either an additional negative-pressure m atter constituent 
withi n General Relativity (Dark Energy), or a modification of General Relativity (Dark Gravity l)Gu and HwaneL 
2002)). This review rests on the observed global homogeneity and isotropy of the universe (Robertson- Walker cos- 
mology RW), and will emphasize the difference between Robertson- Walker kinematics and dynamics. We will recall 
when RW symmetry could determine the cosmodynamics: To explain the observed present cosmological acceleration 
without constituent Dark Energy, Einstein dynamics must be modified at low spacetime (Ricci) curvature. 

Because the homogeneous expansion history H{z) of the global universe measures only kinematic variables, it cannot 
fix the underlying dynamics: cosmographic measurements of the late accelerating universe, are consistent with either 
a static cosmological constant or a dynamic "dark energy" , which itself may be constituent Dark Energy or modified 
gravity (Dark Gravity). The Concordance Model ACDM, with a small static Dark Energy or cosmol ogical constant 
A, gi ves a good fit to all present observations, but also allows a moderately dynamic "dark energy" ijSpergel et all 
2006). We will review how static and dynamic "dark energy" and how dynamic Dark Energy and Dark Gravity could 
be distinguished theoretically and by future cosmological, local or intermediate scale observations. 

Whether Dark Energy or Dark Gravity, "dark energy" has two distinct dynamical effects: it alters the homogeneous 
expansion history H(z) and it suppresses the growth of density fluctuations 5 at large cosmological scale a(i). Because 
the growth function g{z) = S/a depends on both these effects, large angular scale CMB temperature anisotropies, 
the late-time growth of large scale structure, and refined weak lensing observations potentially distinguish static from 
dynamic "dark energy" and Dark Energy from Dark Gravity (Section II). 

In any metric theory of gravitation, the material stress-energy sources T^ v determines the spacetime (Ricci) cur- 
vature tensor or Einstein curvature tensor G^ v = i? M „ — g^R/2. Before considering dynamical alternatives, we 
recall how the asymptotic spacetime curvature, at vanishing matter density, can constrain any RW dynamics, without 
assuming Einstein gravity: Birkhoff's Theorem is a geometric theorem, holdi ng in any spherically symmetric met- 
ric ge ometry whose vacuum has vanishing spacetime (Ricci) scalar curvature l)C. Callan and Peeble 1 1196.3 iPeeblel 
1980). Applied to any spatially homogeneous RW universe, Birkhoff's Theorem asserts that local Newtonian gravity 
fixes the global dynamics of any matter-dominated unive rse whose vacuum is Ricci-flat ()C. Callan and Peebles! . Tl965t 
iMcCrea and Mima . Il934t iMilnel 1 19341 IWeinbergL Il972|) . The Ricci curvature of the vacuum, or cosmological con- 
stant, distinguishes high- from low -curvature alternatives to E i nstein's original gravity theory (Section III. A ). High- 
curvature modifications (such as ijArkani-Hamad et~a U 119981 iRandall and Sundrunl Il999t iBinutrav et all. l2000f) l 
require sub-millimcter correction s to Newton ' s inverse-square grav ity; low-curvature modifications (such as ACDM, 
Dvali-Gabadadze-Porrati (DGP) ijDvali et all l2000i iDeffavetl 1200 if) theories) can preserve Newtonian gravity locally, 
but must be asymptotically dominated by a cosmological constant. Without Dark Energy, our accelerating universe 
is matter-dominated, and Einstein gravity needs low-curvature modification. 

Section III also emphasizes that contrived (epicyclic) dynamic Dark Energy can explain the present acceleration, 
but still cannot explain the Cosmic Coincidence ("Why so small now?"), without fine tuning or anthropic reasoning. 
This will lead us, in Section V, to consider the Dark Gravity dynamical alternatives to Dark Energy. 

Section IV reviews how Dark Energy or Dark Gravity dynamics determines the adiabatic and the effective sound 
speeds, which govern the growth of fluctuations. To illustrate how different dynamics and effective sound speeds can 
underly the same equation of state, we compare canonical (quintessence) and non-canonical (k-essence) scalar field 
descriptions of Chaplygin gas Dark Energy . 

For a truly static cosmological constant, we revert to Einstein's original definition as an intrinsic geometric classical 
parameter. By disconnecting the cosmological constant from energy-momentum sources, we side-step the mysteries of 
why quantum vacuum fluctuations apparently do not gravitate and why the present matter density is rougly equal to 
the present "dark energy" density. We review present local and cosmological constraints on General Relativity, before 
proceeding to prospective cosmological, solar system and isolated galaxy cluster tests for modified gravity. We empha- 
size the difficulties coming tests of relativistic cosmology face: The next decade may dist inguish static from dynamic 
" dark energy" , but will still not distinguish constituent Dark Energy from Dark Gravity ijlshak et all [2005) . Besides 
cosmological tests, low-curvature modifications of Einstein gravity may yet be tested in th e solar system (anomalous 
orbital precess ion, increasing Astronomical unit) or in any isolated rich cluster of galaxies ijLue and Starkmanl 120031: 
Ilorioll2005albl) (Section V). 

In conclusion, cosmological scale modifications of classical Einstein gravity are less contrived than fine-tuned Dark 
Energy and arise naturally in braneworld cosmology. By making intrinsic curvature the source of cosmological 
acceleration, Dark Gravity avoids an additional epicyclic matter constituent, may unify early and late inflation, 
and may be refuted by laboratory, solar system, or galaxy cluster (Section VI). 
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TABLE I Kinematic observables for any RW geometry, in terms of Hubble expansion rate H = a/ a. 



Description 


Definition 


Hubble horizon 
bulk expansion 

conformal time since big bang 

vtvrtvipr motion Hititflnpe VtapW to rerl^Viift 7 

spacetime curvature 


1/H = 1/aH = dri/dN 
da 3 /a 3 = 3dN = 3Hdt = md V 
r,(z) = f*dt'/a(t')=f z °°dz'/H(z') 
A..( z \ - r f z dz' IH(z'\ — r(nn ~ n( z)) 
R = 6(H + 2H 2 ) = 6H 2 {1 + q) 


acceleration 
" slow-roll" parameter 
overall "equation of state 
cosmological jerk 


a/a = H 2 + H, q = d/aH 2 = 1 - dH^/dt = -d(l/Ti)/drj = 1 - e H 

e„ = dH~ l /dt = -dlnH/dN 
w(z) = din #7(1 + z) 3 /3dln (1 + z) = -1 + 2e ff /3 
a /a = H 3 + 3HH + H, j = 'd/aH 3 = 1 + 3H/H 2 + H/H 3 



II. EXPANSION HISTORY H(Z) IN ROBERTSON-WALKER COSMOLOGIES 

Our universe is apparently homogeneous and isotropic (Robertson- Walker) in the large. These Robertson- Walker 
cosmologies are four-dimensional conformally-flat generalizations of General Relativity, in which the spacetime (Ricci) 
curvature R and the Einstein scalar curvature G = 3(k/a 2 +H 2 ) are determined by the matter density p, according to 
the gravitational field equations. The homogeneous expansion of our flat Robertson- Walker universe is described by 
the kinematic (geometric) observables in Table I, wherein the cosmological scale a(t) = 1/(1 + z) and the number of 
e-folds N = In a, so that dN = —d\n(l + z) = Hdt = Tidrj. Overhead dots denote derivatives with respect to cosmic 
time t, so that the conformal Hubble expansion rate 7i = a, the Hubble expansion rate H = a /a, and the Hubble 
time H^ 1 = d dja/cdz is the derivative of the conformal time dM{z)/c back to redshift z. Sub scripts denote pre sent 
values, so that H = 73 ± 3 km/sec/Mpc, H^ 1 = 13.4 ± 0.6 Gyr, cH^ 1 = 4.11 ± 0.17 Gpc ijSpergel et al.U2006|) . 

By measuring the evolution of the mean curvature of the background, cosmography maps the homogeneously 
expanding universe, without reference to dynamics or sources of curvature. However, we will see in Section III.C that 
the asymptotic Ricci curvature or vacuum Ricci curvature Roc = 4x 2 pde{o- — oo) does constrain the Robertson- 
Walker cosmodynamics, distinguishing high- and low-curvature modifications of Einstein gravity. In Einstein gravity, 

G = x 2 p/Z,x 2 = 87tGjv = l/M 2 Pl 

in terms of Newton's constant Gn and the reduced Planck mass Mp. In Einstein's original field equations, H(t) is 
the only degree of freedom, only the tensor components of the metric g^ v are propagating, and (absent a cosmological 
constant) the asymptotic or empty space scalar curvature Roc = 0. When Einstein gravity is modified, H = dH/dt and 
H = d 2 H/dt 2 , or the cosmological acceleration q(t) and jerk j(t) become additional degrees of freedom, describable 
by scalar or vector gravitational fields. 

Conformal flatness means that light propagates in Robertson- Walker cosmologies as in Minkowski space. This 
directly implies a Hubble expansion in cosmological scale a(t), an expansion history H(z), different cosmological 
distances, and other kinematic quantitities listed in Table I. 

A. Kinematics: Distances to Supernovae, Luminous Red Galaxies, Last Scattering Surface 

The CMB shift and the first baryon acoustic peak are standard rulers measuring proper motion distances 

d M (z)= f dz'/H(z') 
Jo 

back to the last scattering surface at redshift z r = 1089 and to luminous red galaxies typically at redshift z\ — 0.35, by 
observing the CMB shift parameter S = V^mO-ffo^Af (^r) = 1.716 ± 0.062 dSpergel et al||2 00fjl')and the first baryon 
acoustic peak A = y/fl m0 H [d 2 M (zi ) jz\ H (zi )] 1/3 = 0.469 ijFairbairn and Goobarl 120051 lEisenstein et alll2005|) . Cali- 
brated supernovae la are standard candles at low redshift z < 1.7, whose observed flux=absolute luminosity/47T(iL(;z) 2 ; 
meas ures their luminosity distance di,(z) = (1 + z)dia{z) l)Perlmutter et al 1 119991: IRiess et all l2004t lAstier et all 
2005). 

These cosmological distances then map t he evolution histo ry H(z), but only after differentiation with respect to 
redshift. Differentiating the Gold Sample ijRiess et all 12004)) supernova luminosity distances dij{z), Tegmark et al 
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l)Wang and Tegmarkl 12005) obtained Figure 1, a plot of h(z) = H(z)/100. A second differentiation of the observed 
distances gives the overall "equation of state" w(z) = 'y(z) — 1 = dh\H 2 /{\ + z) 3 /3dln(l + z), decreasi ng from 
in the matter-dominated epoch, to —2/3 at present, and apparently tending towards —1 in the far future ((TegmarkL 
l2002t IWang and TegmarkL |2004|) . Determining the overall effective "equation of state" thus requires two numerical 
differentiations of the sparse, noisy primary data on cosmological distances. 

As will be discussed in Section III. A below, the late accelerating expansion is conventionally described by an effective 
mixture of General Relativity ideal fluids, now consisting of pressure-free matter (baryons+CDM) and "dark energy" 
defined by p E E = 3MpH 2 — p m and jeie = —dlnpuE/^dN, so that the overall stiffness and "equation of state" 

7 = -d\\iH 2 /2>dN = 7 m O m + Jde^de, ^de = 1 - fi m 
w = w m Q m + w DE Q. DE = w DE (l - Q m ). 

If Dark Energy exists as a matter constituent in General Relativity, then wde is its "equation of state". Otherwise, 
wde{1 — f^m) = d\nH 2 /(l + z) 3 /3<iln(l + z) measures the Dark Gravity modification to the Einstein-Friedmann 
equation. 

Because f2 m o ~ 1/3, the effective "dark energy equation of state" is now wdeo ~ — 1, so that the "dark energy", 
whether Dark Energy or Dark Gravity, is now static or quasi-static, nearly a cosmological constant A = x 2 pde ~ 
2Hq 2 . Indeed, the recent three-year WMAP data, with narrowed constraints on Q m h 2 = 0-126 ± 0.009, O m = 
0.234 ±0.035, n s = 0.961 ± 0.017, expansion ag e t Q = 13.73^ ; " Gy r, large scale structure l|Tegmark et all 12004^ and 
supernova data, makes wdeo = —0- 926io'o75 l|Spergel et all l2"006j) . consistent with the Concordance Model ACDM 
and severely limiting how dynamic any "dark energy" can now be. 



B. Expansion History Does Not Determine Cosmodynamics 



The expansion history H(z) constrains but does not determine the cosmodynamics. Most simply, in Einstein's field 
equations 

G^ v = x 2 T^/i, x 2 = 8nG N = 1/M 2 , 
whose time-time component is the Einstein-Friedmann equation 

k/a 2 + H 2 = K 2 p/i, 

the dynamics and expansion history cannot distinguish between a cosmological constant (static Dark Gravity) — Ag^ 
added to the left side and a constant vacuum energy density (static Dark Energy) p vac = hg^/ x 2 on the right side. 
This familiar Dark Energy/Dark Gravity degeneracy persists when both are dyna mic. Thus, Table II presents five 
two-paramet er fits to the combined Superno va Legacy Survey ijAstier et all l2005|) , baryon acoustic peak and CMB 
observations ijMaartens and Maierottol l2006'). The first and best fit is to the Concordance Model 

(H/H ) 2 = O m o(l + z) 3 + f^Ao + £Ikq> Slmo + Qao + ^K0 = 1- 
The next two fits are to spatially flat Dark Energy models, in which the Einstein-Friedmann equation is 

(h/h q ) 2 = a m0 (i + zf + q xo (i + zf^+ra, n m0 + n xo = 1, 

and the "equation of state" is wde = const or wo + w a z/(l + z), with p ast-average wr>F. = 0-/N) wqe(N') dN' = 
wde = w — w a zln (1 + z)/(l + z) ijLindeil l2005|) . The last two fits ijAstier et all l2005t iMaartens anci Maierr^fo , 
2006) ar e to the original spatially curved Dvali-Gabadadze-Porrati (DGP) Dark Gravity model ijDvali et all 2000; 
lDeffavedl200ltlBento et all 12006^ 

(H(z)/H ) 2 = [1/2(3 + v/(l/2/3) 2 + r! m0 (l + z) 3 ] 2 + Q K0 (l + z) 2 , O m0 + Q K0 + y/l - Q K0 /f3 = 1, 

and to a generalized flat DGP model 

(H/Ho) 2 = n m0 (l + zf + (1 - Q m o)(H/Ho) a , 

fitted only to SNLS+BAO and assuming the prior fi m0 = 0.27+° ^ llFairbairn and Goobarl l2005t lAlam and SahniL 
|2005j) . The spatial curvature £Iko is set equal to zero in the (|F airbairn £uio^ ^baiTl2^05TTAstl^ , 
2005) models and emerges as practically zero in the (Ma artens and Maierottol l2006j) models. The dynamical Dark 
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FIG. 1 The cosmic expansion history (dimensionless Hubble parameter h(z)=H(z)/100) from the Riess et al. (2004) Gold 
sample (top panel) and from simulated future data (bot tom panel) for the NASA/JDEM mission concepts JEDI (solid points) 
and SNAP (dotted points) ( Wa ng and Tegmarkl . |20 05Tl . The existing Gold sample data and the simulated future JEDI are 
both consistent with AC DM (solid curve) (from (IWang and Tegmarkl 12005^ 1. 



TABLE II Model Fits to the Observed Expansion History: Two parameter fits to the joint SNLS, BAO and CMB shift data 
by the Concordance Model, two different dynamical Dark Energy, and two different Dark Gravity cosmological models. 



Model 


n K 


SlmO 


other parameters fitted 


Concordance Model (ACDM) (Maartens and Maierotto, 2006) 


-0.0050 


0.265 


fi A = 0.740 


flat constant wde (Asticr ct al., 2005) 

flat wde(z) — wo + w a z/(l + z) (Lindcr, 2005) 






0.271 ±0.021 
0.260 


wde = -1.023 ±0.087 
w = -0.78, w a = 0.32 


original Dvali-Gabadadze-Porrati (Maartens and Maierotto, 2006) 
generalized flat DGP (Fairbairn and Goobar, 2Q05J 


-0.0297 



0.260 

n 9 7 +0.06 
u ' z ' -0.04 


f3 = H r c = 1.39 
n - _n 17+0.87 

a — "J- 1 ' -0.63 
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Energy and Dark Gravity expansion historie s are equivalent under the substitution (1 + Z ) 3 ( 1 -+ 1 "- D - B ) <-> (H/Hn) a , with 
an average vwe = -1 + a/2 = -l.QSto^ l|Dvali and Turned . 1200.1) . The Gold SN data l|R.iess et all 120041) would 
have yielded slightly poorer fits than the SNLS data we have used. 

Because the supernovae are observed only at low redshifts and the CMB first acoustic peak and the luminous red 
galaxies at recombination redshifts z r — 1089 and z± = 0.35, other smooth parameterizations could fit the past data 
equally well. In any case, because the c osmological distan ces involve two integrations over w(z), they all smear out 
information on the "equation of state" ijMoar et all 1200 This requires sm ooth parametrization of the "equation 
of state" and binning of the sparse, noisy data ( W ang and Teem arkl 120041) , and justifies using no more tha n two 
parameters for present and next decade observations ijLinder and Hutererl 120051 [Caldwell and Lindeit 2005]). (I 11 
retrospect, because observations constrain the directly observable H 2 (z) and the "dark energy" density better than 
its derivative, the "equation of state" , it might have been better to parameterize the past average wde(z), rather 
than w D e{z) l)Wang and FreeseL l2004j) ). 

The observed evolutionary history is already somewhat better fitted by the static ACDM model 
l|Maartens and Maierottol l200fJl than by a ny dynamical Dark Energy or Dark Gravity model in Table II, and the 
latest WMAP data Wdeo = — 0-926^0 075 ^Snergel et all 120061) . even more severely limits how dynamic any "dark 
energy" can now be. Nevertheless, the uncertainties still allow some late-time evolution of w(z). Our purpose will 
now be to discriminate among these nearly sta tic alternatives by observing the fluctuation growth factor on Hubble 
horizon scales. Figure 2, from l(TegmarkL 120021) . shows the ranges of redshift and conformal length scales over which 
such spacetime fluctuations are likely to be measured cosmologically, in the next few years. 

III. CLASSIFICATION OF ROBERTSON-WALKER COSMOLOGIES BY THEIR VACUUM SYMMETRY 

A. Homogeneous Evolution Conventionally Described by General Relativity Perfect Fluids 

By General Relativity, we understand Einstein's original field equations without cosmological constant, before 
introducing the cosmological constant on the left (geometric) side of his original field equations, changing the field 
equations from Einstein to Einstein-Lemaitre and the asymptotic Ricci curvature from flat to curved. This introduces 
into the classical action a very small energy scale A ~ Hq <g; Mp. By avoiding identifying the cosmological constant 
with any right-side matter stress-energy content, this classical approach distinguishes static from dynamic "dark 
energy" and avoids considering the two cosmological constant problems, why quantum vacuum energies apparently 
do not gravitate and why the prese nt matter density p m p ~ M 2 , A . (Replacing t he Einstein Lagrangian R by R — 2 A 
is equivalent to unimodular gravity l)Buchmuller and Dragorlll988tlUnruhl ll989) in which the Einstein-Hilbert action 
is varied holding Jg = — 1: Instead of appearing in the Hilbert-Einstein action, A then enters as an undetermined 
Lagrange multiplier.) 

Although Robertson- Walker cosmology does not assume General Relativity, its expansion history may be expressed 
in terms of an equivalent perfect fluid defined by p = 2>MpH 2 = p m + pr>E and 

overall equivalent pressure: P/c 2 = -Mp(3H 2 + 2H) = P m + P DE 

overall "equation of state": w = Pj pc? = —(1 + 2H/3H 2 ) = — 1 + 2e#/3 = w m Cl m + wde^de 

overall equivalent fluid stiffness: 7(2) = 1 + w = —2H/3H 2 = 2e#/3 = j m Cl m + Ide^de 

overall equivalent enthalpy density: p + P/c 2 = —2M p H = —dp/'iHdt = p m + P m — dpp>E/3dN. 

So defined, pde is either constituent Dark Energy or a Dark Gravity addition to the Einstein- Friedmann equation 
that is non-linear in H 2 . For ordinary non-relataivistic matter p m ~ a~ 3 ,j m = X,P m = = w m - Thus, ever since 
matter dominated over radiation, the overall "equation of state" w — wde{1 — Omo). 

Our universe apparently evolved from a high-curvature de Sitter (early inflationary) universe P = —p = const. 
Assuming the Weak Energy Condition w > —1, so that no phantom matter or cosmic rip intervenes, it will expand 
monotonically H < 0, towards a different (late inflationary) low-curvature de Sitter universe. The deceleration —q(t) 
decreased from 1, when radiation dominated, and is still decreasing, along with the Hubble expansion rate and Ricci 
curvature. 

Deceleration changed over to acceleration when the Hubble horizon changed from expanding to shrinking with 
conformal time, when the "slow-roll" parameter e# = dH^ 1 jdt = —dlnH/dN fell below 1, and w fell below —1/3. 
The acceleration has already increased to present values qo » 0.52, w ~ —0.74, but because e^o w 0.4, this recent 
inflation is still far from truly slow-rolling. The jerk j(t) increases monotonically from a minimum j m i n = —1/8 
when w = —1/2, towards j — > 1, as the universe asymptotes towards a terminal de Sitter universe, with low constant 
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FIG. 2 Shaded regions show ranges of scale and redshift over which various cosmological observations are likely to measure 
spacetime fluctuations over the next few years. The lower left region, delimited by the dashed line, is the non-linear regime 
where rms density fluctuations exceed unity in the ACDM model (from ( Teg markL 120021) '). 



curvature Roo = 12i?^ ~ 4iJg. Such a de Sitter attractor explains the growth of "dark energy", but not why it now 
approximates the energy density of ordinary matter. 

During three long epochs listed in Table III, the universe is dominated by a single barotropic phase with a constant 
equation of state w and diminishing Ricci scalar curvature R(t), in which the scale a(t) ~ ^2/3(i+i«)^ ^ _ 2/3(1 + w)t 
and the acceleration and jerk are fixed at q = — (1 + 3w) /2, j = 1 + 9w(l + w)/2. When these perfect fluids are mixed 
or when cosmological scalar fields appear, the "equation of state" w(z) and the jerk j(z) = l + 9w(l + w)/2~3dw/2dN 
change, the composite fluid is imperfect and supports entropic perturbations. 
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TABLE III Kinematics and Ricci curvature for barotropic phases with power-law growth a ~ ^ 2 / 3 ( 1 +™) j a ^ exp(Ht). 



w 


a(t) 


H(t) 


q(t) 


j(t) 


R(t) 


Model Universe 


1/3 


t l/2 


l/2t 


-i 


3 





radiation 





t 2/3 


2/3t 


-1/2 


1 


3H 2 


non-relativistic matter (Einstein-de Sitter) 


-1/3 


t 


1/t 








6H 2 


curvature energy (constant conformal expansion^, Milne 


-1 


exp Ht 


const 


1 


1 


12H 2 


empty de Sitter 



B. Dark Energy Requires Fine Tuning to Explain the Cosmic Coincidence 



Dark Energy is reviewed in l|Padmanabhanl 12003) . If it exists, Dark Energy is usually attributed to an additional 
ultra-light scalar field 0, so that p = p(p, (f>) is adiabatic only when the scalar field is frozen or tracks the background. 
Defining X = d^^d^ip, the scalar field is canonical (quintessence) when its kinetic energy is X/2, non-canonical (k- 
essence) when the kinetic energy is non-linear in X. Canonical quintessence is driven by a slow-rolling potential and 
can track the background matter, making dw/dz > 0. k-essence is driven by a non-canonical kinetic energy and can 
be arranged to switch from tracking during radiation dominance over towards a dominating cosmological constant, 
making dw/dz < recently. Of course, this different dynamics gives quintessence and k-essence different clustering 
properties. 

The present "dark energy" density pdeo ~ 2HqM p <gc Mp. While ultra- light dynamic Dark Energy can evolve 
down to this very small value, it does not explain the Cosmic Coincidence, "Why so small now?", why pdeq ~ 2/5 m o, 
without the extreme fine-tuning it was invoked to avoid. Canonical and non-canonical Dark Energy ultimately require 
fine-tuning: quintessence, in order to explain the Cosmic Coincidence; k-essence, in order to initiate the transition 
towards a cosmological constant after radiation dominance ends. 



C. Without Dark Energy, Our Universe Must Be Asymptotically de Sitter 

Birkhoff's Theorem l)C. Callan and Peeblesl HoM iPeeblesL 1198(1) : In any locally isotropic (spherically symmetric) 
system whose vacuum is Ricci-flat (i? /iy =0), the vacuum metric must be Schwarzshild: 

9tt =9rr =l-2G N M{r)/r, 

where M(r) is the mass interior to r. For any small spherical shell in empty space, the Newtonian potential must 
vanish inside, and de crease as 1/r outside. Birkhoff's Theorem is a geometric theorem, which generalizes Newton's 
iron sphere theorem l)C. Callan and Peeblesl Il965l IPeeblesL Il980(l from Newtonian gravity to Einstein gravity or any 
high-curvature modification of Einstein gravity. 

Application to any Robertson-Walker Cosmology l)C. Callan and Peeblesl Il965l IWeinberS . Il972j) : Applied to a 
homogeneous universe with matter density p(a), M(r) = ATrp(a)r 3 /3, Birkhoff's Theorem has remarkable dynamical 
consequences. In a homogeneous expanding universe, a small comoving shell lying at r(t) = A*a(t), encloses a mass 
M(r) = A'^47T/9(a)a 3 /3, and has constant Newtonian energy 

r 2 /2 - G N M(r)/r = \l[d 2 /2 + >c 2 p{a)a 2 /6]. 

Using Birkhoff's Theorem, Milne and McCrae JMcC rea and derived the global Friedmann 

equation 

a 2 — x 2 pa 2 /3 = const, k/a 2 + H 2 = x 2 p/2> 

for any pressure-free universe, without assuming Einstein's field equations. If Dark Energy does not exist and the 
vacuum is Ricci flat, Birkhofff's Theorem would make our presently pressure-free universe decelerate according to this 
Friedmann expansion equation. 

Classification of Dark Gravity Theories: If there is no Dark Energy, the Ricci curvature of the vacuum distin- 
guishes high-curvature from low-curvature modifications of General Relativity: Robertson- Walk er universes whose 
vacua remain Ricci-flat in fou r dimensio ns (e.g. Arkani-Hamed et al ijArkani-Hamad e~ a[|,|l998), Randall-Sundrum 
(Ran dall and Sundrumlll999|l . Binutray (Bi nutrav et all 120001) can only modify Einstein gravity in the ultra-violet. 
Robertson- Walker universes which maintain Einstein gravity locally can only modify Einstein gravity cosmologically 
(e.g. ACDM, self-accelerated DGP). If there is no Dark Energy, our accelerating universe is now dominated by 
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pressure-free matter, so that the Friedmann-Robertson- Walker equation must be modified at low Ricci curvature, by 
introducing a cosmological constant. 

When Einstein gravity is modified: If the universe is asymptotically Ricci-curved, the modified Friedmann equation 

k/a 2 + H 2 = x 2 f(p)/3, 

maintains Einstein gravity at high density f(p) — * p, but crosses over to de Sitter — ► const = Mp A at scales r c 
approaching the Hubble horizon H(z) . About any isolated source of mass M and Schwarzchild radius r$ = 2GnM, 
Einstein gravity remains a good approximation only for distances r < r* up to the Vainstein scale r* = (rsr 2 ) 1 ^ 3 ~ 
{HQrsY^H^ 1 <C H^ 1 at which 2GnM/t* = rs/r* = H 2 {a)r 2 — r 2 /r 2 . These are the promising low-curvature 
modification of Einstein gravity that asymptote to de Sitter, to be discussed in Section V.B. 

IV. DYNAMICS DETERMINES THE GROWTH OF FLUCTUATIONS 

Allowing inhomogeneities breaks translational invariance, leading to Goldstone mode sound waves that lower the 
CMB angular power spectrum at large scales (low multipoles I) and leads to growth of large scale structure. In a 
mixture of cosmological fluids or dynamic scalar fields, the equation of state is generally not adiabatic: fluctuations 
propagate in the conformal Newtonian gauge with an effective sound speed c 2 — P.x/p.x — w — dw/3(l + w)dN), 
generally different from the adiabatic sound speed c 2 = dP/dp = P / p. For canonical scalar field quintessence, c 2 = 1, 
but for non-canonical k-essence, c s can vary rapidly, nearly vanishing near the radiation/matter cross-over, where 
w{z) is changing. 

Because the entropic pressure fluctuations are proportional to (1 + w)(c 2 — c 2 ), they are insensitive to the effective 
sound speed in the quasi-static limit w(z) ~ — 1. This minimizes the differences between static and dynamic "dark 
energy" and between any dynamical Dark Energy and Dark Gravity, making their fluctuation growth factors hard to 
distinguish, in present and in next-generation experiments (Section V.C). 

With the same equation of state and adiabatic sound speed, different dynamics generally leads to different effective 
sound speeds. This equation of state degeneracy is illustrated by the toy Chaplygin Gas, whose adiabatic equation 
of state P — —A/p and sound speed c 2 = —w(z) can be derived from adiabatic fluid dynamics, from a non-canonical 
Born-Infeld scalar field, or from a canonical tracking scalar field with potential V(<fi) = (\/A/2)[cosh3</>+ 1/ cosh 30]. 
If derived from a perfect fluid or from the Born-Infeld Lagrangian C$= — VoVl + x 2 X, with non-canonical p = 
Vo/Vl — x 2 X, P = — VbVl — x 2 X, w = — 1 + x 2 X, the perturbations are adiabatic, and the effective sound speed 
c 2 = —w(z) = c 2 . But, if derived from the canonical scalar field with potential V ((/)), entropic perturbations make the 
effective sounds speed c s = 1. Agreeing in the static limit, these three theories give the same adiabatic sound speed 
c 2 — —w(z). Differing dynamically, they show different effective sound speeds, c 2 — —w(z), — w(z), 1. 

This toy model illustrates how the adiabatic sound speed depends only on the equation of state, but the effective 
sound speed and growth of fluctuations depend on dynamics. All three Chaplygin models can fit the kinematic 
observations of the expansion history, bu t the Born-Infeld model fails dynamically by providing insufficient power in 
the observed large scale mass spectrum l|Amendola et all 120031) . This failure can be remedied by generalizing the 
equation of state to P = —A/p a with a ~ 0, so that this generalized C haplygin gas is nearly indistinguishable from 
a cosmological constant l)Amendola et all 120031: iBento etal l l200alZhil|2004h . 

V. DARK GRAVITY: DYNAMICAL MODIFICATIONS OF GENERAL RELATIVITY 

Because Dark Energy is contrived, requires fine tuning and apparently cannot be tested in the laboratory or solar 
system, we now turn to Dark Gravity as the alternative source of cosmological acceleration. This Dark Gravity 
alternative arises naturally in braneworld theories, naturally incorporates a classical extremely low spacetime intrinsic 
curvature, and may unify "dark energy" and dark matter, and possibly early and late inflation. We will see how Dark 
Gravity, besides being tested cosmologically, can also be tested in the solar system, Galaxy or galaxy clusters. 

A. Present Local and Cosmological Tests of General Relativity 

General Relativity is a rigid metric structure incorporating general covariance (co-ordinate reparametrization in- 
variance), the Equivalence Principle, and the local validity of Newtonian gravity with constant Gm, in the weak field 
and non-relativistic limits. General covariance implies four local matter conservations laws (Bianchi identities). The 
Einstcin-Hilbcrt action is linear in Ricci curvature, so that the Einstein field equations are second order, the two 
tensorial (graviton) degrees of freedom are dynamic, but the scalar and vector g^ v degrees of freedom are constrained 



10 



to be non-propagating. Quantum gravity has always motivated high-curvature (Planck scale) modifications of General 
Relativity. Now, the surprising discovery of the accelerating universe motivates extreme low-curvature (cosmological 
scale) modifications of General Relativity. 

General Relativity differs from Newtonian cosmology only by pressure or relativistic velocity effects, which are 
tested in the solar system, in gravitational lensing of light, in the primordial abundance of light elements, in the 
dynamical age, and in the large angular scale CMB and late-time mass power spectrum. Therefore, in order of linear 
scale, modifications of General Relativity must be sought in: 

• labor atory violations of the Equivalence Principle (Eotvos experiments) and solar system t ests llDamour et al. , 
1990) (lunar ranging, deflection of light, anomal ous orbital precessions of the planets, Moon ll Lue and Starkmari , 
l2003t iLui l2005tlGabadadze and Idesiasll2005|) . secular increase in the Astronomical Unit ijloriol 20053)) 

• galaxy and galaxy cluster number counts (|loriol [2005b) 

• gravitational weak lensing 

• cosmological variation of Newton's Gn and other "constants" 

• the suppression of fluctuation growth on large scales or at late times. 



B. Classification of Modified Gravity Theories by Spacetime Curvature of Their Vacuua 

As already noted, in General Relativity only the tensor degrees of freedom in the metric are propagating and 
the homogeneous RW evolution depends only on H(z). If the Einstcin-Hilbert action is modified, additional scalar 
and vector degrees of freedom will appear, and the evolution will also depend on H, H, which can be represented 
by scalar or vector gravitational fields. The basic distinction between high- and low-curvature modifications of 
General Relativity depends on the spacetime (Ricci) curvature of their vacua. It is simplest to begin by considering 
four-dimensional metrical d eformations of General Relativity, which are often inspired by string-theory or M-theory 
ijDamour and Polvakovl,H994albr) or appear as projections of higher-dimensional theories. 



1. Extra Degrees of Freedom in Four-Dimensional Gravity 

• Scalar-tensor gravity, th e simplest and best-motivated extension of General Relativity ijFuiii and Maedal . l2003t 
iGapozziello et al1l2005[) : In the original Jordon frame, a scalar gravitational field proportional to time- varying 
1/Gn, is linearly coupled to the Ricci scalar R. After a conformal transformation to the Einstein frame, the 
scalar gravitational field is minimally coupled to gravity, non-minimally coupled to matter. In the Einstein frame, 
the gravitational field equations look like Einstein's, but the matter field is coupled to the scalar gravitational 
field as strongly as to the tensor gravitational field, so that test particles do not move along geodesies of the 
Einstein metric. Test particles move along geodesies of the original Jordon metric, so that the Weak Equivalence 
Principle holds. 

Scalar-tensor theories modify Einstein gravity at all scales and must be fine-tuned, in order to satisfy obser- 
vational constraints. Nucleosynthesis and solar system constraints severely restrict an y scalar field compo- 
nent, rendering any Dark Gravity effects on the CMB or H (z) evolution imperceptible l)Bertotti et all [2003: 
ICatena et all 12004 ICapozziello et al.l 120064 - 

• higher-order metric f(R) theories: Stability of the equations of motion allows the Lagran gian to depend only 
on R, and only trivially on other curva ture invar i ants, P = R^R^ , Q = Ra^sR"^ 5 ijOstrogradskil [1850) 
or derivatives of any curvature scalar ijWoodardl |20O6). The se f(R) theories are equivalent to scalar ten- 
sor theories with vanish ing Brans-Dicke parameter u>bd = ijTevssandier and Tourrericl Il983t lOlmd . [2005; 
ICapozziello et al.U2006bfh 

The simplest low-curvature modification l)Carroll et all 12004" 120051: iNoiiri and Odintsovl l200r|l . replacing the 
Einstein Lagrangian density by R— n 4 /R, leads to accelerated expansion at low-curvature R < /s 2 ~ Hq, but has 
negative kinetic energies and is tachyonically unstable. This instability would be tole rable in empty space, but 
would be vastly and unaccept ably amplified inside matter llDolgov and Kawasaki |2003|) . and phenomenologically 
unacceptable outside matter ijSoussa and WoodardLl2004|) . These f(R) theories, like their scalar- tensor gravity 
equivalents, can be fine-tuned to avoid these potential instabilities and satisfy supern ova and solar system 
constraints llWoodardl. 120061 ISoussa and Woodardl 12004 INoiiri and Odintsovl I2004l200l . but not cosmological 
constraints l|Amendola et alJ . l2006(l . 
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• TeVeS (relativistic MOND theory): Adding an additional vector gravitational field, could explain galactic 
rotation curves and the Tully-Fisher relation, without invoking dark matter, and possibly unify dark matter 
and "dark energy" ((Bekensteinl . I2004T) . Because gravitons and matter have different metric couplings, TeVeS 
predicts that gravitons should travel on geodesies different from photon and neutrino geodesies, with hugely 
differ ent arrival times fro m supernova pulses. It also predicts insufficient power in the third CBR aco ustic 
peak llSkordis et a,lll200fift . In any case, now that WMAP data requires dark matter l)Spergel et all feOQfil) . the 
motivation for TeVeS disappears. 



2. Extra Dimensional Modifications of Einstein Gravity 

In extra dimensional braneworld th eories, scalar fields appear naturally as dilatons an d modify Einstein gravity 
at high- curvature, by bra ne warping llRandall and Sundruml . Il999l iBinutrav et all 12000^ . or at low- curvature, by 
brane leakage of gravity l)Dvali et all |2000|K If quantized, these theories encounter serious theoretical problems 
(ghosts, catastrophic ultra-violet instabilities, strong coupling problems). Until these problems can be overcome, 
these theories can only be regarded as effective field theories, incorporating an extremely low infra-red scale at 
low spacetime curvature. This suggests an infra-red/ultra- violet connection, since effective field theories ordinarily 
incorporate ultra-violet paramet ers. 

In the orig inal DGP model l|Dvali et, all 1200(1 iDeffavet et all l2002alrl) . the brane's finite stiffness leads to an 
effective modified Friedmann equation, 

H 2 + k/a 2 - H/r c = H 2 + k/a 2 - HH /(3 = x 2 p/3 
(H(z)/H ) 2 = [1/2/3 + V(l/2/3) 2 + fWl + z) 3 ] 2 + fi*o(l + z) 2 , 1 = fi m0 + Sl K0 + V± - W/3 

on the four-dimensional brane by inducing an additional curvature term H/r c at the cosmological scale (3 = Hor c = 
1.39, r c = /3/Hq = H^ 1 ~ 5.7 Gpc. This modified Friedmann equation interpolates between Einstein's pressure-free 
universe at large redshifts, and the empty de Sitter universe with constant Hubble expansion = l/r c = Hq/(3, 
in the asymptotic future. The universe began its late acceleration at z acc = (20 m o//3 2 ) 1 ^ 3 — 1 ~ 0.58. This is the 
original DGP model fit on the fourth line of Table I, which turns out to be spatially practically flat. 
In flat 3-space, the modified Friedmann equation has the self-accelerating solution 



H = v/x 2 p/3 + (l/2r c ) 2 + 1/2 



H(z)/H = 1/2/3 + v/(l/2/3) 2 + r» m0 (l + z) 3 , 1 = Q m0 + 1/(3, 

Because this self-accelerating solution has a Ricci-curved vacuum, Birkhoff 's Theorem does not apply on the 4D brane, 
and Einstein gravity still holds at the shortest distances. However, about any isolated condensation of Schwarzchild 
radius r$ = 2GnM/c 2 , the self-accelerating metric 



g tt = 1 - r s /2r + ^ r|r/2r|, g rr x = 1 + r s /2r - ^r 2 s r/8rl, r<r±, 
and Einstein gravity already breaks down at the Vainshtein scale l)Dvali et all l2003j) defined by 

r* = (tst*) 1 '* - {Hvrsf^H, 1 « H Q \ 

This scale, surprisingly intermediate between rs and Hq , is also where the growth of fluctuations changes from 
Einstein gravity to linearized DGP or sc alar-tensor Brans-Dicke growth, with an effective Newton's constant slowly 
decreasing by no more than a factor two ijLue et al. .120041). 



The original flat DGP model can be generalized I Dvali and Turneil l2003|) to 



H 2 -H a r«- 2 = k 2 P /3, l/!3={l-n m0 ) x l^- a \ l = n m0 +/3 a ~ 2 , 

which is equivalent to a "dark energy" p^E = 3MpH 2 — p = 3MpH a r"~ 2 , wde = — 1 + ct/2. This generalization 
reduces to the original flat DGP form for a = 1, but otherwise interpolates between the ACDM model for a — 0, (3 — 
1.18 and the Einstein-de Sitter model for a = 2, j3 — oo. For small a, it describes a slowly varying cosmological 
constant. Excluding the CMB shift data, the joint S NLS-BAO data fits this nearl y static Dark Gravity model fit with 
a = — 0.17±g;||, (3 = 1.16 or r c = (3/H - 8 Gpc l|Fairbairn and Goobail Eool . This is the generalized flat DGP 
model on the last line of Table I. 
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FIG. 3 The proper motion or conformal distance r(z) — c J* dz' /H(z') back to redshift z, calculated for flat ACDM (dashed) 
and DGP Dark Gravity (solid) models with present matter fraction fl m o = 0.3. The DGP Dark Gravity model is also mimicked 
by a Dark Energy model with w(a) = —0.78 + 0.32z/(l + z). Between static and dynamic "dark energy", the differences in 
distan ces and in expansion history H(z) = cdz/r(z) are small, but the differences in growth factor will be larger in FIG. 4. 
(from iKov ama and Maartcns. I2005T) ) . 
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FIG. 4 The linear growth history g(a) = S/a for flat ACDM (long dashed), DGP Dark Gravity (thick solid) and Dark Energy 
models in Fig. 3. Because the DGP model Newton's "constant" weakens with time, it shows a little more growth suppression 
than that in the mimicking Dark Energy model. DGP-4D (thin solid) shows the incorrect result obtained by neglecting 
perturbations of the DGP 5D Weyl tensor. The 5D Weyl tensor perturbations make dynamical Dark Gravity (DGP) and Dark 
Ener gy hard to resolve, but both dynamical models are distinguishable from static Dark Energy, (from (Kovama and Maartens, 
2005)) 
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C. Prospective Tests of Modified Gravity 

Prospective cosmoloqical ob servations: Figure 3 compares the recent expansion histories Koyama and Maartens 
ama and Maartena l2005[) calculated in the static flat ACDM model with the dynamic flat DGP Dark Gravity 
model (3 = 1.41, chosen to fit Q m o = 0.3, and with its Dark Energy mimic. The degeneracy in these three models 
is lifted by their different linear growth factors plotted in Figure 4: DGP Dark Gravity always suppresses growth a 
little more than Dark Energy does, but substantially more than the Concordance Model. The present normalized 
linear growth factors g(a) — 5/a — 0.61, 0.68.0.80 for DGP Dark Gravity, Dark Energy, ACDM respectively will lead 
to slightly different large-scale CMB and gravitatio nal weak lensing (cosmic shea r) convergence effects. 

These linear growth factors have been calcu lated ijKovama and Maartensl 12005) on subhorizon scales and agree with 
( Ishaketal 1 120051: iLue and Starkmanll2003|) . but are not yet reliable on super horizon scales, where gravity leakage is 
most important. Unfortunately, the large-scale CMB angular power spectrum is obscured by cosmic variance and by 
foreground effects of nearby structures. Nevertheless, the linear growth factors already suggest that next-generation 
observations may distinguish static from dynamic " dark energy" , but will be unable to disting uish Dark Energy f rom 
DGP Dark Gravity. To do so will require a newer, more ambitious weak lensing shear survey (|lshak et all 12005). 

Prospective solar system and galaxy cluster observations: The modifications to Einstein gravity at the Vainstein 
intermediate scale r* may also be tested in next-generation sola r syste m measurements of anomalous precessions of 
plane tary or lunar orbits ijLue and Starkmanl l2003t iDvali et all [2003) or of a secular increase in the Astronomical 
Unit (Iorio, 2005b). These Vainstein scale modifications may also be observable in precision tests about other isolate d 
stars (r$ ~ 3 km, r+ ~ 280 pc) or about isolated spherical galaxy clusters r$ ~ 100 pc, r+ ~ 28 Mpc) florid l2005al) . 

VI. CONCLUSIONS: ACDM, FINE-TUNED DARK ENERGY, OR MODIFIED GRAVITY 

We have reviewed present and prospective observations of "dark energy" as the source of the observed late cosmo- 
logical acceleration, in order to emphasize the differences between kinematical and dynamical observations, between 
static and dynamic "dark energy" , and between Dark Energy and Dark Gravity. We conclude that 

• Cosmological acceleration is explicable by either a small fine-tuned cosmological constant or by "dark energy" , 
which is now at most moderately dynamic. This "dark energy", if dynamic, is either an additional, ultra-light 
matter constituent within General Relativity, or a low-curvature modification of Einstein's field equations. 

• The best and simplest fit to the expansion history, static "dark energy" (ACDM), interprets the cosmological 
constant as a classical intrinsic geometric Ricci curvature, rather than as vacuum energy. This side-steps the 
cosmological constant problems, why quantum vacuum fluctuations apparently do not gravitate and why the 
current matter density is roughly that of "dark energy" . 

• The homogeneous expansion history can also be fitted by moderately dynamical " dark energy" . Only observing 
the large-scale inhomogeneity growth rate will distinguish between dynamic and static " dark energy" . 

• Projected cosmological observations of the growth factor in the large-scale angular power spectrum, mass power 
spectrum, or gravitational weak lensing convergence should distinguish static from dynamic " dark energy" , but 
not Dark Energy from Dark Gravity 

Originally invoked to explain late cosmological acceleration within General Relativity, quintessence and k-essence 
Dark Energy fail to explain the Cosmological Coincidence "Why Dark Energy appears now?" , without fine-tuning or 
anthropic reasoning. Low- curvature modifications of Einstein gravity, such as DGP, are conceptually less contrived 
than finely-tuned Dark Energy, and arise naturally in braneworld theories. They explain cosmological acceleration 
as a natural consequence of geometry, may unify early and late inflation,, and may even be tested by refined solar 
system or galaxy observations. 
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